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Abstract 

Multivariate probability density functions of returns are constructed in order to 
model the empirical behavior of returns in a financial time series. They describe the 
well-established deviations from the Gaussian random walk, such as an approximate 
scaling and heavy tails of the return distributions, long-ranged volatility-volatility 
correlations (volatility clustering) and return- volatility correlations (leverage effect). 
The model is tested successfully to fit joint distributions of the 100+ years of daily 
price returns of the Dow Jones 30 Industrial Average. 
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The methods developed in studying complex physical systems have been suc- 
cessfully applied throughout decades to analyze financial data [1,2,3] and they 
continue to attract gradual interest [4,5,6,7,8,9,10]. The field of research con- 
nected to modeling financial markets and development of statistically based 
real-time decision systems has recently been named Econophysics. In this pa- 
per, we construct a phenomenological model for description of the multivariate 
distribution of returns in a financial time series. 

The random walk model proposed by Bachelier in the year 1900 [1] is equiv- 
alent to the Gaussian multivariate probability density function (PDF) of the 
returns £$: 

The absence of the correlations, Corr[&,£j] = 0, from a time frame greater 
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than ti — ti-i = 20 Min [6] has been widely documented and is often cited 
as a support for the efficient market hypothesis [8]. The multivariate PDFs 
constructed in this work are extensions of the Gaussian PDFs, aimed to model 
the well-established deviations in the behavior of financial time series from the 
Gaussian random walk. 

The Levy stable truncated univariate PDFs [2,9] are known to provide, for a 
financial time series, (i) an approximate scaling invariance of the univariate 
PDFs with a slow convergence to the Gaussian behavior and (ii) the existence 
of heavy tails. We propose the multivariate Student PDFs, 

for modeling the empirical PDFs with returns The marginal PDF (2) is 
again PDF (2). If we integrate out all of the & except for one, we get (2) with 
n — 1. The tails of the distributions behave empirically like [6] ~ d£/£ 4 , and 
so a ~ 3. 

For the Student PDF (2), we have E[&] = and Corr[|&| , |&|] = 2/(2 + tt) = 
0.39. The correlation of the returns vanishes as in the Gaussian random walk. 
The square of the volatility equals of = Var[£j] = uofa/(a — 2). 

The cumulative returns, £ = are described by 

Tl 71 > 1 > 

dw(o= I dm-E^s:^)Ud^ = suh4 ( 3 ) 
i=i 1=1 

where Q 2 = X^Li oof. The variance of the £ increases linearly with n, in agree- 
ment with the Gaussian random walk and in the rough agreement with the 
empirical observations. Eq.(3) represents the scaling law for financial time 
series. 

The multivariate Student PDFs have therefore heavy tails and the exact scal- 
ing invariance from the start. These distributions can be modified further to 
describe two other well-established stylized facts which are (iii) long ranged 
volatility- volatility correlations that are also known as volatility clustering [11] 
and (iv) return-volatility correlations that are also known as leverage effect 
[12,13]. 

The empirical facts show that there is a slow decay of the correlation func- 
tion. An extension of the PDF (2) that has the value Corr[|£j| , which is 
decaying with time is rather straightforward. We use uncorrelated multivariate 
distributions for different groups of the returns. The analogy with the Ising 
model can be useful: The groups & with the same multivariate Student PDF 
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can be treated as domains of spins aligned in the same direction. We assign 
the usual probability for every such configuration 



n-l 



w[<Ji 
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where <7j = ±1. The normalization constant is given by 1/N = 2(2 cosh(/9)) n_1 . 
The correlation of the absolute values of the returns equals 2/(2 + n) provided 
that £j and belong to the same domain. Otherwise the result is zero. The 
probability of getting the £j and within the same domain can be found to 
be Wi = e~ l1 where e -7 = e _/3 / [eP + e"' 3 ) < 1 and I = i — k. The coefficient 
Corr[|£j| , for the modified PDF has therefore the form 2/(2 + tt)wi. The 
absence of the correlations would formally correspond to (3 = +oo (7 = +00). 
This is the case when the multivariate PDF is a product of the n univariate 



In order to incorporate leverage effect, we consider the cu-function depending 
on the signs ej_ p = sign(^_ p ) of the lagged returns (p = 1,2, ...). The val- 
ues €i = ±1 are assumed to be independent variables which take the two 
values ±1 with the equal probabilities. The q dependence is modeled by 
uji = C (l — Y^=\ uPe i-p) where C is a normalization factor. Negative 
recent returns €i- p = — 1 (p = 1,2, ...) increase the volatility, so p > 0. The 
value of p is connected to the overall strength of leverage effect. 

Taking volatility clustering and leverage effect into account, we obtain for 
I = i - k > 0: 



where g\ = E[uiUk} and hi = E[ui€kUJk] can be calculated analytically. 

The empirical correlation functions for the 100+ years of daily price returns 
of the Dow Jones 30 Industrial Average are fitted using a superposition wi = 
Cl e-^ 1 + c 2 e" 72 ' for I > with c x = 0.18, c 2 = 0.08, 71 = 1/1200, and 
72 = 1/233 (else c = 0.74 and 70 = +00, so that X^ =0 c m = 1). We have also 
used p = land v = exp(— 1/16). The results are in a good agreement with 
the data. The quantitative comparison can be found in Ref . [14] . The value of 
c i + C2 = 0.26 indicates that 26% of the empirical PDF consists of products 
of the multivariate Student PDFs. The value Corr[^, |£ fc |] vanishes at i > k in 
agreement with the observations, since the uj{ depend on lagged returns only. 
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Corr[|&| , \i k \\ 
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The complete multivariate PDF of the returns is given by 



2 n-l 



exp(-/3 m (n - 2s - 1)) 
(2cosh(/5 m ))"- 1 



S%(£ n , ■■;£l)c— E C ™ E 



s+1 




(7) 



where = n + 1, n > n s > ... > ri\ > 2, and n = 1. The marginal 
probability of the PDF (7) is a PDF (7) again. 

The authors wish to thank the Dow Jones Global Indexes for providing histori- 
cal prices of the Dow Jones Averages. M .I.K. is grateful to Metronome- Ricerca 
sui Mercati Finanziari for kind hospitality. 
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♦ Empirical deviations from the Gaussian random walk: 

(i) approximate scaling, 

(ii) heavy tails of the return univariate distributions, 

(iii) long ranged volatility-volatility correlations (volatility clustering), 

(iv) return- volatility correlations (leverage effect). 

♦ Multivariate Student PDFs for modeling (i) - (iv). 

♦ Conditional and marginal PDFs. 

♦ Pricing Log contracts. 
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Why Econophysics? 



Methods developed in studying complex physical systems have 
been successfully applied through decades to analysis of financial 
data 

L. Bachelier, Ann. Sci. Ercole Norm. Suppl. 3, 21 (1900), 
P. Levy, Theorie de /'Addition des Variables Aleatoires 

(Gauthier-Villars, Paris, 1937), 
B. B. Mandelbrot, J. Business 36, 294 (1963), 
R. N. Mantegna and H. E. Stanley, Nature 376, 46 (1995) 
and continue to attract gradual interest. 

Last years: 

Development of statistically based real-time decision systems 
results in the growing importance of the Econophysics research. 

Over 10% of the papers in Physica A are related to 
Econophysics. 

New journals: Quantitative Finance, ... 
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Generalized Wiener process 



Financial time series are described by a generalized Wiener 
process: 

dS(t) 



S(t) 



= /idt + df(t). 



Here, d£(t) is a noise with the expectation value and variance 

E[d£(t)] = 0, Var[d£(t)] = <r(t) 2 dt. 

The volatility a(t) is a measure of the magnitude of market 
fluctuations. 

We consider a discrete random walk and set 

dt = ti - ti_i, Si = S(ti), & = d^(ti), di = cr(ti). 
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^ Gaussian random walk 

The random walk model proposed by Bachelier in 1900 year is 
equivalent to the Gaussian multivariate PDF: 



No statistically significant correlations at ti — U-i > 20 Min 



The multivariate PDFs constructed in this work are extensions of 
the Gaussian PDFs to model the well established deviations in the 
behavior of financial time series from the Gaussian random walk: 

(i) heavy tails of the return distributions, 

(ii) approximate scaling, 

(iii) volatility clustering, and 

(iv) leverage effect. 



G„(6) = , 1 /0 exp(--Va i - 2 A 2 ) ff - 
VS,y (27r) n / 2 FV 2^ 1 i = \(Ti 
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Deviations from Gaussian random walk: Autocorrelations 



Autocorrelation function of £(t) and autocorrelation function of 
\£(t)\, with sampling time interval At = 1 Min. 
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Results are from: Yanhui Liu, Parameswaran Gopikrishnan, Pierre Cizeau, 
Martin Meyer, Chung-Kang Peng, H. Eugene Stanley, 

E-print arXiv: cond - mat/9903369 (1999). 
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Deviations from Gaussian random walk: Heavy tails 

Probability density of At = 5 minutes increments of S&P 500 
index future prices compared to the Gaussian distribution for the 
same mean and variance. The data set contains 75 000 prices. 



Probability density of price changes 

S&P 500 Futures, 1991-95 
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Results are from: R. Cont, E-print arXiv: cond — mat/ '97 '0507 '5 (1997). 
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Deviations from Gaussian random walk: Heavy tails 



Student PDF 
Gauss PDF 



Comparison with the empirical data: 
100+ years of daily quotes of the DOW Jones 30 Industrial 

Average 

The a = 3 Student and 
Gaussian univariate PDFs 
as compared to the hys- 
togram for the DJIA. 

The value of £ is the in- 
crement (noice) added to 
the path followed by the 
index value, a 2 = Var[£], 
The common scale is 
fixed by fitting the vari- 
ance. The solid line de- 
scribes the Student PDF, 
the short-dashed line cor- 
responds to the Gaussian 
PDF. 
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Deviations from Gaussian random walk: Leverage effect 



The leverage effect in the DOW Jones 30 Industrial Average from 

1900 until 2000 year 
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Results are from: Josep Perello, Jaume Masoliver, 

E-print arXiv: cond - mat/0111334 (2001) 
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Student random walk 




Levy stable truncated PDFs proposed by Mantegna and Stanley 
(1995) are known to provide (i) approximate scaling and (ii) heavy 
tails. Two other stylized facts (iii) volatility clustering and (iv) 
leverage effect . Therefore 

MULTIVARIATE STUDENT PDF 

to model realistic distribution of returns: 



Q;+n 
2 



(a7r) n / 2 r(f 



a-\-n 
2 




The marginal Student PDF is again a Student PDF. If we 
integrate out all the increments ^ except for a one, we get (#) 
with n = 1. The tails of the distributions behave empirically like 

~ <i£/£ 4 , and so a ~ 3. 
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Student random walk: Heavy tails 



We compare: 

♦ The n = 1, a = 3 Student PDF, 

♦ The Gaussian PDF, and 

♦ The empirical univariate PDF. 

NB: The empirical PDF is constructed for the 100+ years of the 
daily quotes of the DOW Jones 30 Industrial Average for a period 
of 26.05.1896 - 31.12.1999 (totally 28507 working days) 

The drift rate and variance are found to be 



The PDFs are shown as functions of the £/<r with a 2 = Var[£]. 



H = 0.000257 & Var[£] = 0.000117. 



November 15, 2003 

APFA4, Warsaw, Poland 



Return Distributions in Financial Time Series (page 10) 
with application to pricing derivatives and risk management 



E. Alessio, V. Frappietro, 
M. I. Krivoruchenko, L. J. Streckert 



Metronome-Ricerca sui Mercati Finanziari, 



Torino 



r 



Deviations from Gaussian random walk: Heavy tails 



Student PDF 
Gauss PDF 



Comparison with the empirical data: 
100+ years of daily quotes of the DOW Jones 30 Industrial 

Average 

The a = 3 Student and 
Gaussian univariate PDFs 
as compared to the hys- 
togram for the DJIA. 

The value of £ is the in- 
crement (noice) added to 
the path followed by the 
index value, a 2 = Var[£], 
The common scale is 
fixed by fitting the vari- 
ance. The solid line de- 
scribes the Student PDF, 
the short-dashed line cor- 
responds to the Gaussian 
PDF. 
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Student random walk: Autocorrelations 



For Student PDF, we have 

E[&]=0 & Var[&] = uo] 



a 



a-T 

while the autocorrelation functions have the form 



Corr 


= o, 


Corr[|&| 


= o, 


Corr[|6|,|6|] 


2 

2 + tt' 



The square of the volatility equals 



2 2 GL 
(J: = LU: 

1 a-2 



The correlation function of |^| is positive and does not depend on 

the lag I = i — j. 
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Student random walk: Exact scaling 



If all components of the vectors i/j 
~ 7V(0, 1), a random vector 



i, ...,^ n ) and 77 = (771, ...,7^) 



6 = ^i^i 



with 77 = |r/| has the Student PDF. 
The cumulative effect, £ = YJl=\ii, is described by 



dw(o = / wit - Eeos^te) n d 6 = s ?(^ 

i=l i=l 



ft 



where = £)™ =1 o;f . The value £ can be represented as 



£ = (cJi^i + 002^2 + ... + CJ n ^n) 



2 rsj 



flip 



a 



7]' 



where ip ~ 7V(0, 1). The variance of the £ increases linearly with n, 
in agreement with the empirical observations. We got the scaling 
law, observed empirically by Mandelbrot (1963). 
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Student random walk: Conditional PDFs 



The Student conditional PDF at n^> I = n — k that gives a 
forecast density / steps ahead: 



where 



1 fc 



(i) 



The conditional PDF of the values £ n , is close to the 

Gaussian PDF. 



The ARCH models: The increments are distributed as ~ 7V(0,i; 2 ) 
with u being a function of the lagged increments. The estimator 
(1) is one of the possible estimators quantifying the volatility 
within the ARCH framework. 
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Student random walk: A summary 



The multivariate Student PDFs have heavy tails, the ex- 
act scaling invariance, and volatility clustering from the 
start. They are apparently an excellent starting approxi- 
mation for precise modeling the empirical PDFs. 



These distributions can be modified further to describe two other 
well established stylized facts which are 

(iii) decaying long ranged volatility-volatility correlations (volatility 
clustering) and 

(iv) return-volatility correlations (leverage effect). 
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Decay of correlations 



The long ranged correlations occur, since the value rj in the 
representation 



t i a - 



is common for all the increments Decay of the correlations 
take different 77's for different groups of the 

ANALOGY WITH ISING MODEL: 

The components ^ with the same denominator rj 2 can be treated 
as domains of spins aligned in the same direction. The usual 
probability: 





n-l 


w[(7i, ...,cr n (3 


'] = Nexp(-/3 ^i^i+i) 




i=i 



where cr,- = ±1. 
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Decay of correlations 



The normalization constant: 



l/N = 2(2cosh(/3)) 



n-l 



The correlation of |^| and |^| equals the old value provided that 
these increments belong to the same domain and zero otherwise. 

The probability to get the ^ and within the same domain: 




where 



e _|_ e -0 



and 1 = 1 — k. Corr[|^| , |^|] has the old form multiplied by wi 
Superposition with different 7 is possible. 
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Leverage effect 



To incorporate the leverage effect =^> uji depending on the signs 

€i_ p = sign(&_ p ) of the lagged increments (p = 1,2, ...). The 
function uji = x(e«_i,e«_ 2 , ...) > is defined by 



= C 1 - 



1 - v 00 

P 22 ^ P£ i-P 

" P=l 



We require < z/ < 1 to ensure the convergence of the series. 

Negative recent returns 6^_ p = — 1 (p = 1,2, ...) increase the 

volatility, so p > 0. 

The normalization condition: E[^] = 1, so C = (1 + p 2 



The value of p gives the overall strength of the leverage effect, 

while log(l/^) is the decay time. 
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Results 



Taking the volatility clustering and the leverage effect into 
account, we obtain: 

Cc»r[| 6 | ,|&|] = g a-SWjfo-l) 

7T 9o-^2 



Corr[|&| ,£ k 



2(l->/ + f 



where = E^cj^], /i/ = Ef^e^^], and I = i — k > 0. We have 



21 *% = (! + / ^ 2 



1 + 1/ 



1 + 1/ 



) 2 + 4p 



(2) 



(3) 



?!-!/) ,4 1/ 2 (1 - I/) 3 2/ 



l + v 



(1 + 



o 1 — Z/ N , \ — V i 

{l + p 2 —- )h = -2p i/ 1 



l + i/ 



v 
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Deviations from Gaussian random walk: Autocorrelations 



The empirical correlation functions for the daily quotes of the DJIA are fitted 



Wl = c x e~ lxl + c 2 e~^ 1 I > 



with a = 0.18, c 2 = 0.08, 71 
70 = +oc, so that J2m=0 °m 



0.15 
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0.00 



-0.05 



-0.10 




: 1/1200, and 72 = 1/233 (else c = 0.74 and 
1), also p = 1 and v = exp(— 1/16). 

Correlation functions ver- 
sus time lag I = i — k be- 
tween two increments & 
and £fc. The empirical 
data are constructed us- 
ing the daily quotes of the 
DOW Jones 30 Industrial 
Average. The volatility 
clustering (circles) and the 
leverage effect (boxes) are 
clearly seen. The data 
are compared to predic- 
tions of the Student multi- 
variate PDF shown by the 
solid lines. 
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More about Student random walk 



REALISTIC MULTIVARIATE STUDENT PDF 



^ _ ^ „ ^ exp(-/3 m (n - 2s - 1)) 

m=0 s=0 



(2cosh(/3 m )) n - 1 



x 



s+l 

^ II ^n f -n f _i (^n f -l 5 
n s ,...,n f=l 



£n f _i) 




where = n + 1, n > n s > ... > n\ > 2, and n = 1. The 
calculation time for computers increase exponentialy with n. The 
straightforward calculation of the n = 13 PDF takes ~ 1 Min. The 
marginal probability of the PDF (9) is a PDF (9) again: 



n 



, ... , 



£i)f II d£i = S£(&,...,£i)F 

i=k+l 
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More about approximate scaling 



The value of c = 0.74 turned out to be unexpectedly large. It 
indicates that 



Three quarters of the empirical PDF consist of the product of 

the Student univariate PDFs. 



The product does converge to the Gaussian PDF due to the 
central limit theorem. One should thus verify that the proposed 
PDF nevertheless supports the approximate scaling observed for 
financial time series. 





A qualitative explanation: the Levy stable PDF fits well the 
central part of the empirical distributions and therefore of the 
Student PDF, providing an exceptionally slow convergence to the 
Gaussian PDF. 
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More about approximate scaling 



The realistic PDF is constructed for the daily increments. There 
are no a priori reasons to believe that the daily scale is more 
appropriate than other time scales. 



The criterion for selection of the best time interval is 
the quality of the fit of the empirical univariate PDF. 



It is seen from Fig. 4 that e.g. At = 13 days would be an 
inaccurate choice. At the same time, the excellent fit of the daily 
empirical PDF indicates that the assumption 



At = 1 day 



on which the realistic PDF is based is quite reasonable. 
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Deviations from Gaussian random walk: Approximate scaling 

Scaled univariate Student PDFs 



10 u 

6 
5 
4 
3 



10 1 

6 
5 
4 

3 



2 



10 



2 - 



10" 



Student PDFs 
Gauss PDF 




The a = 3 scaled Student 
PDF for n = 1,2,3,5,8, 
and 13 (solid lines) as com- 
pared to the Gaussian PDF 
(short-dashed line). The 
lower vales of the n cor- 
respond to the higher val- 
ues of the PDF at the ori- 
gin £ = and to the lower 
values at the high absolute 
values of the £. The con- 
vergence to the Gaussian 
PDF with increasing the n 
is slow. The value of £ 
is the noise added to the 
path followed by the index 
and a 2 is variance of the 
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Student PDF n = 13 
Student PDF n = 1 
Gauss PDF 



Deviations from Gaussian random walk: Approximate scaling 

Scaled n = 13 distributions based on the realistic modified 
multivariate Student PDF as compared to the empirical PDF for 

100+ years of the DJIA daily quotes 

The n = 13 and n = 1 Stu- 
dent and Gaussian PDFs as 
compared to the hystogram for 
the 100+ years of quotes of 
the DOW Jones 30 Industrial 
Average. The value of £ is 
the increment (noice) added to 
the path followed by the index 
value, a 2 = Var[£]. The com- 
mon scale of the distributions 
is fixed by fitting the variance 
of the £. The solid line stands 
for the n = 13 modified multi- 
variate Student PDF, the long- 
dashed line stands for the n = 1 
Student PDF, and the short- 
dashed line denotes the Gaus- 
sian PDF. 
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Conditional vs maginal PDFs 



Ensemble of uncorrelated stocks marginal probabilities 



Single index =>* time averages => conditional probabilities 



The above discussion refers to the marginal probabilities. If no 
correlations exist, the marginal probabilities coincide with the 
conditional probabilities. 



Ergodic hypothesis is valid, since all the val- 
ues 7*. and \w(\jv) are distinct from zero, so 
the correlation lengths are all finite. 



It means that the time averages are equivalent to the ensemble 
averages, provided that the time scale used for the estimation is 
much longer than the largest correlation length (I/71 ~ 5 calendar 
years). 

At smaller time scales, however, the distinction between these two 
^Jbypes of the probabilities can be significant. 
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Conditional PDFs 



The conditional PDF of the increments £ n , £1 equals 

S„Kn, £i)f 



Sn(Cn, .-j £k+l|£k 3 £i)f — 



Sk(^k, £i)f 



It is normalized to unity. For i > k 

Var[^|^,...,^i] F 



(5) 



o-(i-fe) v exp(— /3 m (z — i; — ^)) 

2. 2^2, (2c OS h(/3 m ))-«+- 



e iv? e fc+i m=0 



(6) 



where ft = min(l, v — 1), and Ki = max(0, fc + 1 — u). 
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Conditional PDFs 




Var[^i|^i_i, ...,£i] F 



(7) 



while the old expression 



can be referred to as the marginal (unconditional) volatility. The 
ensemble average of (12) gives (13). 

It is hard to distinguish empirically between the increments £ ?: and 
returns A^. The asymptotic dW ~ d^/^ 4 ~ e _3A G?A implies 
E[S n ] < oc,and Var[S n ] = oo, while dW ~ dA/A 4 ~ ^/(£ln 4 £) would 
cause more severe theoretical problems: E[£ n ] = VarfS^] = oo. The 
asymptotic behavior ~ dt;/^ 4 refers naturaly to the increments 
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Increments vs returns 



Return on an asset per time interval [t^, 



A^lnSi/Si-! 



-2 



nd the drift rate //are small, and so 



The increments ^ ~ 10 

~ & to a first approximation 

The Ito's lemma can be applied to lnS(t) provided that 
E[(d£(£)) n ] < oo for all n. In our case, the asymptotic distributions 
of the increments are such that E[(o?^(t)) 2n ] = oo for n > 2. 
Equation 



Ai = ln(l +fi + &) 



cannot be simplified further by expanding the logarithm in power 
series of the 

v - / 
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Log contracts 



The multivariate PDF of returns: 



n 



dW[Ai] = s«(6) F nd6 

i=l 



where ^. 



The knowledge of the multivariate PDF provides, in prin- 
ciple, the most complete information on the future market 

behavior. 



A simple check: Fair value of a market index at t = tk 



n 



E[S n ] = S k E[exp( £ A i )|£ k ,...,£ 1 ] F = (l + r) 1 S k 

i=k+l 

where ji = r is the risk-free discount rate corresponding to the 
expiration date t n . The result has the conventional form. The 
dispersion of S n , however, diverges: Var[5 n |£ fc , = 00 
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Log contracts 



Pricing of the Log contracts: The fair value of a long position, 

L = lnSn/Sfc, is defined by: 



n 



E[L] = E[ ]T Ai|a,...,Ci]F 

i=k+l 

Using the conditional PDF, we obtain: 



n 



n 



E[ £ A i |e k ,...,? 1 ] F ^2- 1 £ £(r--a? k ). 

i=k+l e n ,...,e k+ i i=k+l 

The result depends on the lagged increments through 

crf k = Var^il^k, ...,fi] F . 

The pricing of the Log contracts is well defined both with respect 
to the fair value and its dispersion. 
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Conclusion 



CONCLUSION: 

The asymptotic behavior is important for the correct estimates 
of financial risks. The use of the multivariate Student PDF 
gives more precise idea on the involved risks and, in particular, 
on those ones connected to large market fluctuations. 

The following four stylized facts 

( 



i) heavy tails, 

ii) approximate scaling of the returns, 

iii) volatility clustering, 

iv) leverage effect 

are well established empirically, but missing in the Gaussian 
random walk model. We constructed the multivariate PDFs of 
incremenets and returns for financial markets, that take 
approximately those four effects into account. 

♦ The model can be useful for more accurate pricing derivative 
securities and risk management. 
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